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INTRODUCTION 


In  the  previous  reporting  period'  we  reported  on  measurements  of 
the  optical  distortion  that  was  observed  when  a CVD  zinc  selenide  window 
was  exposed  to  a 2 kWatt  C.W.  laser  beam  at  10.6  pm.  A very  preliminary 
theoretical  expression  was  given  to  account  for  the  observed  time -dependent 
distortion  when  a window  is  simultaneously  heated  by  a laser  beam  and  cooled 
by  gas  flow  across  the  faces. 

It  was  recognized  that  a more  exact  theory  would  be  useful  not  only 
for  predicting  OPD  but  also  for  evaluating  the  level  of  stress  development. 
However  an  exact  theory  of  the  time  dependent  effects  for  an  arbitrary  laser 
beam  profile  has  not  been  developed  in  closed  form.  The  method  most 
frequently  used  is  to  model  the  window  and  the  heating  rate  by  a mesh  of 
lumped  parameters  and  then  apply  nodal  point  analysis  by  a computer  to 
analyze  the  temperature  and  stress  profiles. 

This  is  a satisfactory  method  for  specific  experiments,  but  it  is  also 
dcsireable  to  have  an  analytical  solution  so  that  the  relations  between  the 
individual  system  parameters  may  be  more  easily  understood  and  the 
scaling  laws  perceived.  Fortunately,  the  C.  W.  laser  beam  available 
for  the  experiment  is  very  nearly  gaussian  in  shape.  It  was  found 

possible  to  derive  an  analytical  expression  for  the  time  dependent  temperature 
and  stress  profiles  for  an  infinite  plane  window  exposed  to  a gaussian  laser 
beam  and  simultaneously  to  face  cooling.  This  work  and  the  comparison  with 
an  experiment  was  presented  at  the  ARFA  Laser  Window  Conference  at  Boulder, 
Colorada  on  July  12,  1976.  The  same  material  is  also  given  in  this  report 
and  is  the  principle  work  done  under  this  program  during  the  current  reporting 
period. 


Zar,  J.  L. , Investigations  into  the  Feasibility  of  High  Power  Laser  Window 
Materials,  Semiannual  Technical  Report  No.  3,  April  1976.  To  be  published. 

"Zar,  J.  L. , Semiannual  Technical  Report  No.  2,  Investigations  into  the 
Feasibility  of  High  Power  Laser  Window  Materials,  AFCRL-TR-75-0605, 
Oct.  1975. 


The  transmission  of  a laser  beam  through  a plane  window  gives  rise 
to  a change  in  optical  path  length  as  a result  of  three  physical  phenomena; 

(1)  change  in  window  thickness  by  heating  and  thermal  expansion;  (2)  change 

t 

in  index  of  refraction  as  a result  of  heating  and  (3)  change  in  both  thickness 
and  index  as  a result  of  thermoelastic  stress  set  up  by  temperature  gradients 
during  the  time  of  illumination  or  subsequent  cooling  of  the  window.  In  this 
paper,  a solution  to  the  optical  path  difference  (OPD)  resulting  from  these 
three  mechanisms  is  derived  in  analytic  form  and  the  results  applied  to  a 
window  experiment  where  the  time  dependent  OPD  for  a plane  window  was 
experimentally  determined  during  and  after  illumination  by  a laser  beam. 

The  problem  of  the  time  dependent  heating  of  a laser  window  by  an 

isotropic  and  centrally  positioned  laser  beam  has  already  been  solved  by 
(123)* 

various  methods.  ’ ’ In  the  present  work  the  Fourier  Bessel  trans- 

form is  used  and  it  yields  the  temperature,  strain,  and  stress  distributions 
in  closed  form.  The  window  response  can  thus  be  evaluated  without  the 
need  for  numerical  integration  methods,  at  least  for  some  problems  of 
practical  importance. 

A typical  laser  window  material  is  isotropic  with  failure  characteri- 
stics generally  described  as  brittle.  For  an  acceptable  window,  the  absorption 
of  laser  radiation  in  it  is  sufficiently  low  that  the  transmitted  beam  may  be 
approximated  by  an  intensity  distribution  independent  of  z.  The  window 
has  a radius  b,  a thickness,  £.  The  incident  laser  beam  that  is  used  is 
normal  to  the  window  surface  and  has  a gaussian  intensity  distribution 
given  by 


I 


,,  2 2 
t-,  c -a  r / 
P a e /7T 


11) 


The  first  two  authors  use  finite  Fourier  Bessel  transforms  to  state  the 
problem  nut  use  numerical  methods  for  the  integrations.  The  third  author 
does  not  employ  a gaussian  laser  beam. 
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where  P = laser  power  and  a is  a constant.  The  window  has  an  absorption 
coefficient  of  (3  per  cm  ( |3  t « 1)  so  that  the  heating  rate  in  degrees  per 
cubic  cm  per  second  is  given  by 


= T e 
o 


2 2 

■a  r 


and 


pPa 


(Z) 


(31 


o irpc 

where  pc  is  the  specific  heat  per  unit  volume  of  the  laser  window  (assumed 

to  be  a constant).  Thus  far  nothing  has  been  said  regarding  reflections 

from  the  incident  and  emerging  surfaces  of  the  window.  This  should  be 

considered  because  (3  is  defined  in  terms  of  the  intensity  of  radiation 

inside  a medium.  If  one  calculates  the  reduction  in  intensity  resulting 

from  surface  reflections,  and  also  the  increase  in  absorption  resulting 

from  multiple  internal  reflections  of  the  transmitted  radiation,  the  two 

(4) 

effects  compensate  each  other.  One  finds  that  the  absorbed  intensity 
A =.  Iq  pf  independent  of  the  reflection  coefficient  and  the  reflected  intensity 
R % 2 RjIq  independent  of  (3  (R^  is  the  reflection  coefficient  from  a single 
optical  surface).  Thus,  nature  has  cooperated  in  this  instance  in  making 
it  unnecessary  to  formally  include  the  micro  details  of  the  reflection  and 
transmission  process.  All  that  is  required  in  Eq.  (3)  is  that  Pf  be  small 
and  that  the  two  surfaces  of  the  window  be  not  accurately  perpendicular 
to  the  wave  front  in  order  that  the  localized  interference  fringes  resulting 
from  internal  reflections  are  spaced  much  closer  than  l/a  apart. 

The  equation  leading  to  the  temperature  distribution  in  the  window, 
as  described,  is  the  Fourier  Heat  Conduction  Law 

.2, 


ft  T 


* ± 

r or 


_1_  ft_T 
K ft  t 


+ ^ 
K 


- 0 


(4) 


0 ^ r < b 

0 < z < fi 


where  k is  the  thermal  diffusivity  and  the  last  term  is  the  heat  generated 

^ T 

per  unit  volume.  The  boundary  conditions  are  that  -x — =0  at  z-  0 and  f; 

3 T ^ z 

and  -jj-p-  = 0 at  r = 0 and  b.  The  case  where  surface  cooling  occurs  will 

be  treated  subsequently  as  an  extension  of  the  solution  of  Eq.  (4). 
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The  differential  equation  is  solved  by  two  calculations  in  which  the 
temperature  at  the  center  is  first  a step  function  and  then  the  solution  is 
rewritten  for  a continuous  central  temperature.  The  first  step  makes 
use  of  the  Hankel  transform  of  order  zero  (also  called  the  Fouriti*  Bessel 
transform)  T (£)-/”  r JQ  (r£)  T (r,t)  dr. 

T(4)  = A f r J (r£)  exp  (-a2r2)  dr  (5) 

° */o  0 

where  Aq  is  the  temperature  at  the  center  of  the  window,  r = 0.  The  inte- 
gration is  given  in  closed  form  by^ 

f(|)  = (Ao/2a2)  exp  (-e2/4a2) 

After  integration,  the  solution  is  found  by  taking  the  inverse  transform 


T(r,  t)  = f e"^  * T(^)  i Jo  (^r)  d| 

Jo 

(6) 

= ^AQ/2a2^  J exp  ^-|2  £l  + 4 k a2tj  /4a2  ^ | JQ  (£  r)  d| 

(7) 

T(r.t) 

= ^Aq/J\  + 4 k a2tj^  exP  f-r^a2/^!  + 4 k a2tj^ 

(8) 

The  second  calculation  uses  the  Duhamel  superposition 
theorem' ^ 

i- 

T(rrt)  = j Ao(t)~  T(r(t-  T)dr 

(<>) 

and  since 

A (t)  = / T dr  = T fdT  = T (t)  . 
° Jo  ° ° 

T(r.t)  = T f (t-JjH  dr 

° Jo  1+4  k a (t-  T) 

(7) 

This  integrates  to 

(10) 

T,„„  - (V4«  .*)[«(- A*) 

Ul) 
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Here  Ei(x)  is  the  exponential  integral  of  negative  arguments,  tabulated 
in  Jahnke  and  Emde*8'  as  the  function  Ei(x)  = - f*  t"  1 exp  (-t)  dt 
(Abramowitz^  calls  this  function  -E^fx)}.  *'X 

Equation  (11)  gives  in  closed  form,  the  time  and  radial  dependence 
of  the  heating  of  a window  by  a laser  beam  of  gaussian  intensity  profile 
when  the  window  is  larger  than  the  laser  beam,  there  is  negligible  heat 
loss  from  surfaces  and  edges  and  there  is  negligible  reduction  in  laser 
intensity  as  a result  of  internal  absorption.  Figure  1 is  a plot  of  the  non- 
dimensional  temperature  rise  T/Tq  t against  the  non-dimensional  window 
radius  ra,  as  given  in  £q.  (11). 

Window  Stress  and  Strain 

The  stress  and  strain  in  a window  transmitting  a laser  beam  of 

gaussian  intensity  can  also  be  stated  in  closed  form.  The  basic  equations 

(9) 

are  those  for  solid  cylinders.  If  the  window  is  thicker  than  the  laser 
beam  mean  radius,  one  should  use  the  long  cylinder  expressions  (plane 
strain)  rather  than  the  thin  plate  formulae  (plane  stress)  because  the  latter 
theory  ignores  axial  stress.  Since  failure  of  brittle  materials  occurs  where 
the  shear  is  a maximum,  the  thin  plate  theory  leads  to  the  conclusion  that 
windows  are  stronger  than  they,  in  fact,  are  because  the  axial  stress  terms 
are  the  largest  of  the  principle  stresses.  The  equations  for  plane  strain 


M - v ) fl  /ixF, 


b 

'^L 


Trdr  - ( 


1 /rZ)f  Trt 

Jo 


f 1 -v  ) o „,/q  E - ( l /b  ) f Trdr  + ( l/r 


’rdr  - T 


( 1 -v  ) ct  /a  E 
z 


( 2v  A>  ) / Trdr  - T 


In  the  above,  a is  the  linear  coefficient  of  expansion,  E is  the 
elastic  modulus,  v is  Poisson1  s ratio;  a is  strain  and  2b  is  the  diameter 
of  the  window.  These  expressions  contain  yTrdr  which  fortunately  is 
integrable  from  Eq.  (11). 
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Figure  4 Non-Dimensional  Axial  Stress 
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From  which  the  individual  stress  values  may  be  obtained  for  an  experi- 
ment by 

S ..  PPciE 


r,  9,  z 8rr»cpc  f 1 - v ) 

It  can  be  seen  from  these  curves  that  the  axial  stress  is  about  twice  as 
large  as  the  others. 

Tho  axial  strain  for  a thin  cylinder  (plane  stress)  is  given  by 
€z  = a T i aj E - v (<rr  + o0)/E 
The  plane  stress  equations  are 


(21) 


(?  2) 


Or  = aE  j^T,/2  - (1/r2)  ^ Trdrj 
o0  = aE  j^/2  + ^l/r2^  f Trdr  - tJ 


a = 0. 
z 


(23) 


These  lead  to  the  expression  for  axial  strain  when  the  laser  beam  nearly 

completely  fills  the  window 

€ = aT  + va  (T  - T.) 

z 1 

For  cylinders  that  are  thicker  than  the  laser  beam  radius  one  uses  the 
stress  Eqs.  (12)  to  (14)  but  in  addition,  one  should  rewrite  Eq.  (22)  by 
replacing  E by  E/(  1 -v  2);  v by  v ./(1-v  ) and  a by  a ( 1+v  ),  This  leads 
to  the  thick  window  e.xpression  for  axial  strain. 

“ -rif  [T  + » 'T  - T,']  <*» 

In  terms  of  non-dimensional  quantities,  one  has  the  non-dimensional  axial 
strain 


4r(l  -v  ) pc k 
' (ltv  )2  ^ * 


V 

v +1 


1 1 T 

a2b2  T t 


(4»ca2t) 


(25) 


In  early  stages  of  window  heating,  the  second  term  in  the  bracket  may  be 
neglected  compared  to  1 so  that 

£ « " — (4/ca^t)  f 2 6> 

Z T t 
o 

This  is  plotted  in  Figure  5 for  various  values  of  the  non-dimensional  radius 
and  transmission  time, 

Failure 

In  triaxial  loading  of  brittle  materials,  failure  occurs  when  the 

maximum  shear  stress  equals  half  the  yield  stress,  a , as  determined  in 

, . . . . (10)  y 

a uniaxial  test. 

There  are  three  maxima  of  shear  which  occur  at  45  degrees  to  the 
principle  stress  directions.  We  expect  failure  when 

T = ± | Max  (Oj  - a2)  = \ oy 

where  o.  and  o_  are  any  of  the  three  principle  stresses  taken  in  whatever 

i b 

order  gives  the  largest  r.  Since  the  axial  stress  is  approximately  twice 

the  other  stresses,  failure  is  predicted  when  a s*  2a  and  this  first  occurs 

z y 

on  the  central  axis  of  the  laser  beam.  The  max  shear  stress  is  plotted  in 
Figure  6. 

Optical  path  changes  in  a plane  window  that  is  heated  by  a laser 

beam  was  observed  by  including  the  window  in  one  leg  of  a Mach  Zehnder 
(4) 

interferometer.  The  dependence  of  the  OPD  on  the  temperature  is 
found  by  evaluating  the  individual  contributions  of  several  phenomena. 

T . T 

OPD  = An!  = f ng  -|~r  dT  + J £g  dT  127) 

where  ng  is  the  index  of  refraction  as  customarily  measured  at  zero 

stress,  £ is  the  window  thickness  and  £ is  the  thickness  measured  at 

s 

zero  stress.  Through  the  use  of  Eq.  (24), 


d£  = « d£z 
dT  ~ dT 


1 - V 


a 

s 


fZH) 
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Here,  a is  the  coefficient  of  expansion  as  measured  at  zero  stress.  It  is 
to  be  noted  that  Eq,  (28)  already  includes  the  thermal  expansion  term  even 
in  the  absence  of  temperature  induced  strains.  Also, 


dn 

dT 


8 n 
da 
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dT 
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9a, 


da 
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9 n 
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da 
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dT 


dn 

s 

dT 


(29) 


0 ~ z 

The  last  term  in  Eq.  (29)  is  required  because  there  is  a change  in  index 
of  refraction  with  temperature  even  when  the  temperature  is  uniform 
and  there  is  zero  thermoelastic  stress.  Following  the  usual  classical 


derivation  of  stress  birefringence,  one  has  for  cubic  crystals 
mitting  a plane  wave  in  the  z direction. 
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= qllaz  + q12(°r  + V 
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(30) 


Here  q 


11 


and  q^  are  the  stress  optic  coefficients. 
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The  derivatives  of  the  stresses  are  obtained  from  Eqs.  (17)  to  (19). 
The  evaluation  of  these  stress-optic  terms  is  very  simple  for  cubic 
crystals  because  the  B term  is  not  present  in  either  a or  in  the  sum 
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Note  that  the  last  term  in  the  bracket  of  Eq.  (32)  is  often  the 
dominant  term. 
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Surface  cooling  of  windows  through  the  use  of  a flushing  gas  stream 
is  an  effective  method  of  removing  laser  beam  generated  heat.  This  is 
particularly  useful  for  windows  with  antireflection  coatings  because  a large 
part  of  the  window  heating  may  originate  in  the  surface  and  it  is  efficient 
to  remove  it  at  the  source. 

Equation  (4)  is  now  expanded  to  include  a term  for  removal  of  heat 
from  two  surfaces  with  a surface  heat  transfer  coefficient,  h W/cm  deg. 
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The  solution  to  this  equation  proceeds  in  a similar  way  to  that  for  Eq.  (4): 

that  is,  by  first  finding  a solution  for  a steady  step  addition  of  heat  by  a 

gaussian  source,  (T/k  = 0 in  Eq.  (33))  and  then  applying  the  Duhamel 

superposition  integral  to  calculate  the  result  for  a time  dependent,  con- 

(8) 

tinuous  heat  source.  ' Equation  (33)  becomes 


where  h'  = 2h/k£, 
to  Eq.  (8)  is 
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The  solution  for  the  step  heat  addition,  corresponding 


TQ(r,  t)  = £t’  /(l  + 4 Ka2t)j  exp  [-fch't]  exp  £-r2a2/(l  4 4 fca2t)j  (35) 

• 'c 

In  the  above  T is  a temperature  change  at  r = 0 which  will  be  evaluated 
shortly.  This  equation  is  now  inserted  in  Eq.  (9) 


“ 

T(r.t)  = i’7  «?■£  dr  (36) 

Jo  1+4  KaT  (t-T) 

The  integration  is  performed  by  expanding  the  first  exponential  in  a power 

series  expf-«h'  (t-T)}  « l-Kh'  (t-T)  + ..  . . For  the  experiment  that  we 

will  try  to  explain,  the  Kh't  « 1 for  t « 100  sec  and  therefore  only  the  first 

two  terms  of  the  expansion  are  sufficient  at  early  times.  This  gives  two 

terms  to  the  integral  of  Eq.  (36)  T(r,  t)  = T^  + T-,.  Tj  is  identical  to 

Eq.  (10)  for  which  the  solution  is  stated  in  Eq.  (11).  This  also  permits 

identification  of  T as  identical  to  T . 

o 

T^  can  be  found  in  closed  form  by  integration  by  pa^ts. 
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T^Ei  [-r2a2J  -Ei  [-r2a2/(l  4-  4 Ka2t)jj  | 

The  complete  solution  of  Eq.  06)  is  thus  obtained  by  adding  Eq.  (11) 


° (4Ka2)‘ 


+ (14  r2a2‘ 


T(  r,  t)  " 7 

4Ka 


j[i+hWA1[Ei(.rV)  - Ei(-^4-Y 

(L  4a2  JL  \ 1 +4fca  t /_ 


4a**  L 


2 2 2 
exp  ( - r a ) - ( 1 + 4 fca  t)  exp 


( 37) 


where  K h1  t " 1.  On  the  axis,  (r  - 0)  this  equation  must  be  restated  to 
remove  the  indeterminate  form. 
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r=0  4fca 
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4a 

Comparison  with  an  experiment.  In  an  experiment,  using  a CW  laser 


beam  at  the  Laser  Window  Test  Facility  a window  of  uncoated  CVD 
zinc  selenide  was  exposed  for  10  seconds  to  a 2 kW  laser  beam  with  an 
approximate  gaussian  intensity  distribution.  A detailed  description  of 
the  experimental  procedure  is  being  published  elsewhere.  In  brief, 

the  distortion  of  the  window  produced  by  heating  was  observed  by  taking 
50  frame  per  second  motion  pictures  of  Mach- 7’ehnder  interference  fringes 
during  and  after  exposure  to  the  laser  beam  using  a helium  neon  laser 
source  to  illuminate  the  interferometer.  During  and  subsequent  to  the 
exposure,  the  two  faces  of  the  window  were  surface  cooled  by  flat  jets 
of  nitrogen  gas  with  an  estimated  surface  heat  transfer  coefficient  of 
0.  0124  W/ctt/  deg  C.  Other  pertinent  data  pertaining  to  this  experiment 
are  given  in  Table  I. 

l'.y  analysis  of  the  motion  pictures,  a curve  of  fringe  count  versus 
exposure  time  was  obtained  for  the  OPD  on  the  center  line  of  the  laser 
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TABLE  I 


EXPERIMENTAL  DATA  FOR  THE  CVD  Zinc  Se  TEST 


Window  diameter 

Window  thickness 

Laser  beam  power 

Beam  mean  diameter 

Surface  heat  transfer  coefficient 

He -neon  laser  wavelength 

Nominal  absorption  coefficient 
(12) 

Poisson's  ratio 

(12) 

Elastic  modulus' 

(12) 

Index  of  refraction 

Temp.  coef.  of  n^1^ 

(12) 

Stress  birefringence 

Volumetric  specific  heat^1^ 
Thermal  diffusivity^ 1 ^ ^ 

( 13) 

Thermal  expansion  coef. 

(13) 

Yield  stress' 


2b 

= 

3.  6 cm 

£ 

= 

. 82  cm 

P 

= 

2 kW 

1/a 

= 

. 24  cm 

h 

= 

.0124  W/cm^  deg 

A 

= 

6.33x10  ^ cm 

P 

= 

.005  cm"1 

V 

= 

.28 

E 

s 

6.  9 x 1011  d/cm^ 

n 

2.403 

s 

dn  /dT 
s/ 

- 

7.9  x 10“5  deg"1 

qll 

= 

-1.39x10  1 ^ cm^/dyne 

ql  2 

= 

0.  58  x 10  1 ^ cm^/dyne 

pc 

= 

2.  65  j/cm1  deg 

K 

= 

. 0491  cm^/ sec 

a 

= 

7.  2 x 10"6  deg"1 

cr 

- 

5.  6 x I011  dyne/cm^ 

y 
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beam.  This  is  plotted  in  Figure  7.  Superimposed  on  this  curve  is  a plot 
of  the  theoretical  Eq.  (32)  evaluated  by  use  of  the  data  in  Table  I.  The 
window  was  supplied  to  us  by  AFCRL  and  it  was  stated  that  the  previous 
measured  value  of  f3  was  0.  005/cm 

The  correction  term  to  the  laser  power,  for  surface  cooling, 
according  to  Eq.  (38)  is  evaluated  first.  After  9 seconds,  it  is  0.985 
and  at  shorter  times  it  is  even  closer  to  1.  0.  Therefore,  it  is  considered 
to  be  negligible. 

The  OPD  is  found  from  Eq.  (32).  The  first  term  in  the  bracket 
(thermal  expansion)  is  1.  18  x 10  deg  V The  second  term  (stress  optic 
effect)  is  13.  9 x 10  ^ deg  *.  The  third  term  (change  in  index)  is 
5.  69  x 10  ■*  deg  * and  is  the  dominant  one. 

Discussion.  Comparison  of  the  calculated  and  experimentally  mea- 
sured curves  of  Figure  7 shows  that  the  theory  accurately  predicts  the 
OPD  up  to  the  first  5 seconds  of  heating.  After  10  seconds  the  theoretical 
curve  is  about  15%  below  the  observed  curve  and  has  a smaller  slope. 

There  are  three  possible  explanations  for  this  discrepancy.  The  central 
portion  of  the  window  has  been  heated  about  25  degrees  in  10  seconds. 

It  is  possible  that  some  of  the  properties  of  CVD  zinc  selenide,  besides 
ng,  have  a temperature  dependence  which  is  not  accounted  for  in  the  cal- 
culations. The  most  likely  property  to  have  a temperature  dependence 
is  p,  which  in  the  present  theory  is  assumed  to  be  constant.  Skolnik, 
et  al.  , calculated  the  temperature  dependence  of  p on  the  assumption  that 
it  is  due  to  multiphonon  absorption  and  stated  that  there  could  be  a 33% 
increase  in  (3  over  a 25  degree  temperature  change.  However  the  authors 
measurements  of  (3  for  the  samples  which  they  possessed  was  independent 
of  temperature  over  the  range  from  300  to  500  °K.  Therefore,  it  is  not 
likely  that  a p dependence  on  T is  the  correct  explanation. 

Another  possible  cause  of  the  discrepancy  at  long  heating  times 
is  that  the  present  theory  assumes  an  infinitely  large  window  with  a cen- 
trally heated  region.  In  actuality,  the  laser  beam  center  was  about  0.  7 cm 
from  the  edge  of  the  window.  The  heated  zone  would  have  diffused  to  the 
edge  in  5 seconds  after  which  time  the  window  temperature  would  rise 
faster  in  this  experiment  than  for  a much  larger  window.  The  extent  of 
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Comparison  of  Experimental  and  Calculated  Optical 
Path  Difference 
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the  extra  heating  may  be  estimated  by  interpolation  of  Figure  l at  ra  = 0 
2 

and  4k  a t = 22  and  44  corresponding  to  5 and  10  seconds  of  heating.  The 
temperature  of  a large  window  should  have  risen  by  50%  for  a 100%  increase 
in  the  time.  For  a smaller  window,  with  a smaller  diffusion  area,  after 
5 seconds  there  would  be  a smaller  increase  in  temperature^  This  is 
qualitatively  sufficient  to  explain  the  difference  between  the  observed  and 
theoretical  heating  curves  at  long  heating  times.  The  third  explanation 
is  that  the  value  of  a may  not  have  been  accurately  measured.  We  have 
calculated  that  for  the  conditions  of  the  experiment,  dT/T  = -3.  24  da  /a  . 
Thus  a 15%  error  in  T could  be  accounted  for  by  a 5%  error  in  the  know- 
ledge of  a.  We  believe  that  a portion  of  the  discrepancy  at  long  heating 
times  may  be  ascribed  to  this  cause  also. 
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FUTURE  PLANS 


f I 

U 


We  have  been  asked  to  assist  Dr.  J.  L.  Stanford  of  the  Naval 
"Weapons  Center  by  providing  exposures  of  some  window  blanks  of  ZnSe, 

KC0  and  NaCff  . Before  doing  this,  we  will  expose  some  ZnSe  plates  to 
the  C.W.  laser  beam  and  compare  the  observed  OPD  with  the  predictions 
of  the  theory. 

It  is  desireable  that  a mechanical  shutter  mechanism  be  used  to  switch 
the  laser  beam  into  the  window,  rather  than  relying  on  the  laser  electrical 
power  and  timing  controls.  A limer  operated  shutter  has  been  installed  and 
it  will  be  used.  It  is  also  desireable  that  the  laser  beam  profile  be  reihccketl 
to  determine  whether  it  is  still  gaussian  and  whether  it  is  circular.  This 
is  also  planned  as  a preliminary  to  the  NWC  window  tests. 
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